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Abstract. We compute the emission of gravitational radiation from pulsating white dwarfs. This is done by 
i using an up-to-date stellar evolutionary code coupled with a state-of-the-art pulsational code. The emission of 

. gravitational waves is computed for a standard 0.6 Mq white dwarf with a liquid carbon-oxygen core and a 

hydrogen-rich envelope, for a massive DA white dwarf with a partially crystallized core for which various 1 = 2 
\f~} , modes have been observed (BPM 37093) and for PG 1159-035, the prototype of the GW Vir class of variable 

stars, for which several quadrupole modes have been observed as well. We find that these stars do not radiate 
sizeable amounts of gravitational waves through their observed pulsation g-modes, in line with previous studies. 
We also explore the possibility of detecting gravitational waves radiated by the /-mode and the p-modes. We find 
that in this case the gravitational wave signal is very large and, hence, the modes decay very rapidly. We also 
discuss the possible implications of our calculations for the detection of gravitational waves from pulsating white 
+3 ■ dwarfs within the framework of future space-borne interferometers like LISA. 
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$_j ■ 1. Introduction 2000) , gravitational wave astronomy will probably be soon 

I a tangible reality. 
Gravitational waves are a direct consequence of General 

Relativity. Many efforts have been done so far to detect Despite of its potential interest, the emission of grav- 
them, but due to the intrinsic experimental difficulties in- itational waves by pulsating white dwarfs has been In- 
volved in the process of detection and in the data anal- tie explored up to now. In fact, apart from the pioner- 
ysis no definite result has yet been obtained. Supernova ing work of Osaki & Hansen (1973), only the gravita- 
core collapse, binary systems involving compact objects, tional wave radiation of rotating white dwarfs undergo- 
and pulsating neutron stars are, amongst others, promis- ing quasi-radial oscillations has been studied so far — see 
ing sources of gravitational waves — see Schutz (1999) for Benacquista et al. (2003) and references therein. White 
a comprehensive review on the subject. Moreover, with dwarfs are the most common end-point of the evolution 
the advent of the current generation of terrestrial grav- of low- and intermediate-mass stars. Hence, white dwarfs 
itational wave detectors, like LIGO (Abramovici et al. constitute, by far, the most numerous stellar remnants 
1992), VIRGO (Acernese et al. 2004), GEO600 (Willke in our Galaxy, outnumbering neutron stars. Moreover, 
et al. 2004), or TAMA (Takahashi et al. 2004), and of the relative simplicity of their physics allows us to ob- 
space-borne interferometers like LISA (Bender et al. 1998, tain very detailed models which can be ultimately com- 
pared with their observed properties. Among white dwarfs 
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* Member of the Carrera del Investigador Cientffico y as ZZ Ceti (or DAV, with hydrogen-rich envelopes and 
Tecnologico, CONICET, Argentina. T off ~ 12 000 K), V777 Her (or DBV, with helium-rich en- 
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velopes and T c g ~ 25 000 K) and GW Vir stars (or vari- 
able PG 1159 objects, with envelopes which are rich in car- 
bon, oxygen and helium, and T c g ranging from ~ 80 000 
to 150 000 K), which show periodic variations in their light 
curves — see Gautschy & Saio (1995; 1996) for reviews. 
The typical periods are within ~ 100 s and ~ 2 000 s and, 
consequently, lay in the region of frequencies to which 
LISA will be sensitive. The luminosity changes of these 
variable stars have been successfully explained as due to 
nonradial (/-mode pulsations. At present, there is a gen- 
eral consensus that variable white dwarfs are very inter- 
esting targets for pulsational studies. Their very simple 
internal structures allow us to predict theoretically the 
pulsational frequencies with a very high degree of detail 
and sophistication. Also, they have a very rich spectrum 
of frequencies which may give us information about the 
stellar mass, the core composition, the mass of the sur- 
face helium and hydrogen layers (if present), the angular 
speed of rotation and the strength of the magnetic field — 
see, for instance, Pfeiffer et al. (1996) and Bradley (1998, 
2001), amongst others. Consequently, it is not surprising 
that in the last years ZZ Ceti and V777 Her white dwarfs, 
as well as GW Vir stars, have been the preferred targets 
for the network called "Whole Earth Telescope" (WET). 
WET observations have been of an unprecedent quality, 
and in some cases have allowed us to disentangle the in- 
ternal structure and evolutionary status of several white 
dwarf stars by applying the powerful tools of asteroseis- 
mology (Nather 1995; Kawaler 1998). 

BPM 37093 is the most massive pulsating white dwarf 
ever found (Kanaan et al. 1992). It is a massive ZZ Ceti 
star — that is, with a hydrogen-rich atmosphere — with 
a stellar mass of ~ 1.05 Mq, and an effective tempera- 
ture T cff ~ 11800 K. BPM 37093 has been thoroughly 
studied (both theoretically and observationally) because 
presumably it should have a sizeable crystallized core 
(Winget et al. 1997). Hence, for BPM 37093 we have de- 
tailed models (Montgomery & Winget 1999; Corsico et 
al. 2005) and extensive observational data (Kanaan et 
al. 2005). Interestingly enough, one of the most appar- 
ent modes of BPM 37093 has a period P = 531.1 s, very 
close to frequency of maximum sensitivity of LISA and 
pulsates with £ = 2. It should be noted at this point that 
£ — 1 modes do not radiate gravitational waves and that 
£ = 2 modes are relevant for the emission of gravitational 
waves, thus making BPM 37093 a especially suitable tar- 
get for LISA. Moreover, the distance to BPM 37093 is 
known (d = 16.8 pc). Consequently, a detailed study of 
the possibilty of detecting the gravitational waves emit- 
ted by this star is of the maximum interest, but still re- 
mains to be done. On the other hand, PG 1159-035, the 
prototype of the GW Vir class of objects, has a complex 
spectrum with several £ = 2 modes (Winget et al. 1991). 
Unfortunately there is not any reliable parallax determina- 
tion for PG 1159-035. More specifically, Werner, Heber & 
Hunger (1991) provide d ~ SOO^qq pc, whereas Kawaler 
& Bradley (1994) obtained d ~ 400 ± 40 pc. However, a 
spectroscopic determination of its mass (M* ~ 0.54 Mq) 



is available. These are, to the best of our knowledge, the 
only two known white dwarf pulsators with confirmed 
quadrupole g-modes. 

In this paper we compute the emission of gravitational 
waves from pulsating white dwarfs. We first compute the 
gravitational waves radiated by a typical 0.6 Mq white 
dwarf with a carbon-oxygen core and a 10 M* hydro- 
gen envelope, which we regard as our fiducial model. For 
this model white dwarf we first compute the gravitational 
waves emitted by (/-modes. Then we compute the gravita- 
tional waves emitted by BPM 37093 and PG 1159-035, the 
only two known white dwarfs with quadrupole (/-modes. 
As it will be shown below, we find that the fluxes radi- 
ated away by these two stars in the form of gravitational 
waves are very small. This is the reason why we also ex- 
plore other possibilities. In particular we also compute the 
fluxes radiated by the /- and p-modes, independently of 
the lack of observational evidence for these modes in pul- 
sating white dwarfs. The paper is organized as follows. In 
§2 we briefly review the basic characteristics of nonradial 
pulsation modes. In §3 we discuss the numerical codes used 
to compute the nonradial pulsation modes of the white 
dwarf models presented here. Sect. 4 is devoted to ob- 
tain the expressions for the emission of gravitational waves 
from pulsating white dwarfs. Finally in §5 we present our 
results whereas in §6 we summarize our findings and we 
draw our conclusions. 

2. Nonradial pulsation modes 

While we highly recommend the interested reader to con- 
sult the excellent monographs of Unno et al. (1989) and 
Cox (1980) for details of nonradial stellar pulsations, we 
believe that a brief overview of the basic properties of 
nonradial modes it is worthwhile to be done. Briefly, non- 
radial modes are the most general kind of stellar oscilla- 
tions. There exist two subclasses of nonradial pulsations, 
namely, spheroidal and toroidal modes. Of interest in this 
work are the spheroidal modes, which are further classified 
into g-, f- and p-modes according to the main restoring 
force acting on the oscillations, being gravity for the g- 
and /-modes and pressure for the p-modes. 

For a spherically symmetric star, a linear nonradial 
pulsation mode can be represented as a standing wave of 
the form % Am (r, 9, <j>, t) = % Am (r) Y^A) e"*-'-™*, 
where the symbol " ' " indicates a small Eulerian perturba- 
tion of a given quantity <J> (like the pressure, gravitational 
potential, etc) and Y™(9, <p) arc the corresponding spher- 
ical harmonics. Geometrically, £ is the number of nodal 
lines in the stellar surface and m is the number of such 
nodal lines in longitude. In absence of any physical agent 
able to remove spherical symmetry (like magnetic fields or 
rotation), the eigenfrcquencies <Jk,i,m are dependent upon 
£ but are 2£+l times degenerate in m. Finally, t m (r) 
is the radial part of the eigenfunctions, which for realis- 
tic models necessarily must be computed numerically to- 
gether with (Tk,e,m- The index k (known as the radial order 
of the mode) represents, in the frame of simple stellar mod- 
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els (like those of white dwarf stars which we shall study 
below), the number of nodes in the radial component of 
the eigenfunction. Generally speaking, g-modes are char- 
acterized by low oscillation frequencies (long periods) and 
by displacements of the stellar fluid essentially in the hor- 
izontal direction. At variance, p-modes have high frequen- 
cies (short periods) and are characterized by essentially 
radial displacements of the stellar fluid. Finally, there is a 
single /-mode for a given £ (> 2) value. This mode does 
not have any node in the radial direction (k = 0) and pos- 
sesses an intermediate character between g- and p-modes. 
In fact, its eigcnfrcqucncy lies between that of the low or- 
der g- and p-modes, and generally slowly increases when £ 
increases. For g-modes (p-modes), the larger the value of 
k is the lower (higher) the oscillation frequency. 

3. Numerical codes 

We compute the nonradial pulsation modes of the white 
dwarf models considered in this work with the help of the 
same pulsational code described in detail in Corsico at al. 
(2001a, 2002). The code, which is based on a standard fi- 
nite differences scheme, provides very accurate oscillation 
frequencies and nonradial eigenfunctions, and has been 
employed in numerous works on white dwarf pulsations 
- see, for instance, Corsico et al. (2004) and references 
therein. The code solves the fourth-order set of equations 
governing Newtonian, linear, nonradial stellar pulsations 
in the adiabatic approximation following the dimension- 
less formulation given in Unno et al. (1989). To build up 
the white dwarf models needed for our pulsational code 
we employed the LPCODE evolutionary code described 
in detail in Althaus et al. (2003, 2005). Our evolutionary 
code contains very detailed physical ingredients. A full de- 
scription of these physical ingredients is out of the scope 
of this paper and, consequenlty, the reader is referred to 
Althaus et al. (2003, 2005) for an extensive description of 
them. Instead, we will only summarize here the most im- 
portant inputs. For instance, the equation of state includes 
partial ionization, radiation pressure, ionic contributions, 
partially degenerate electrons and Coulomb interactions. 
For the white dwarf regime, we include an updated version 
of the equation of state of Magni & Mazzitelli (1979). The 
code uses OPAL radiative opacities — including carbon- 
and oxygen-rich compositions — for arbitrary metallic- 
ity from Iglesias & Rogers (1996) and molecular opacities 
from Alexander & Ferguson (1994). High-density conduc- 
tive opacities are taken from Itoh et al. (1994) and the 
references cited there, whereas neutrino emission rates are 
those of Itoh et al. (1996), and references therein. It is im- 
portant to mention at this point that the stellar models for 
BPM 37093 and PG 1159-035 discussed below have been 
derived from full evolutionary calculations that take into 
account the history of the progenitor stars — see Althaus 
et al. (2003, 2005) for details. Finally, it is worth noting as 
well that during the white dwarf cooling phase, the effects 
of time-dependent element diffusion have been considered 
in the calculations. 



4. Gravitational waves from a pulsating white 
dwarf 

Although the basic formalism for deriving the gravita- 
tional wave radiation of pulsating objects (cither white 
dwarfs or neutron stars) is well known — see, for instance, 
Osaki & Hansen (1973) — we consider appropriate to sum- 
marize it here for the sake of completeness. Moreover we 
will extend it to the case in which a white dwarf has a 
partially crystallized core. Generally speaking, the am- 
plitude of a gravitational wave emitted from any slow- 
moving source in the quadrupole approximation is given 
by (Misner et al. 1973) 



h 



TT 



2G - TT 



(1) 



where "TT" stands for the traceless-transverse gauge, d 
for the distance, and Q is the quadrupole moment of the 
mass distribution, which is defined as 



R3 



p(r)(3xiXj 



S ijr 2 )d d r 



(2) 



As previously stated, we assume that the spatial and tem- 
poral behavior of the the perturbed density profile is pro- 
vided by the following expression: 



P (r, t) = Po (r) + p'(r)Re (¥{"(9, ^)e i<rt ) 



(3) 



where p is the unperturbed density profile, p'(r) stands 
for the radial perturbation of the density profile, Y™(6, <fr) 
are the spherical harmonics, and a is the pulsational fre- 
quency. As we are dealing with £ = 2 modes, and since 
the emission of gravitational waves in this case will be the 
same for all the values of m, we shall choose the simplest 
case. That is, we adopt £ = 2 and m — 0. Additionally, 
it must be taken into account that BPM 37093 has 
a sizeable crystallized core. Therefore, the appropriate 
boundary conditions differ from those of an ordinary 
star. Particularly, the boundary condition at the stel- 
lar center (when crystallization has not yet set in) is 
that given by Osaki & Hansen (1973). However, when 
the core of the white dwarf undergoes crystallization we 
switch the fluid internal boundary conditions to the so- 
called "hard sphere" boundary conditions (Montgomery 
& Winget 1999). Within this approximation the nonra- 
dial eigenfunctions are inhibited to propagate in the crys- 
tallized region of the core. Consequently, and keeping in 
mind that for the axisymmetric case Qn = Q22 = — 5Q33 
and Qij = if i 7^ j (Osaki & Hansen 1973), the following 
expression can be easily obtained 



!33 



111 P 7T 



Jo 

A J 



16tt 

r 4 p'(r) cos(ai) dr 



sin 6>(cos 2 6 - 1)(3 cos 2 6 - l)d6 



(4) 
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Fig. 1. Run of the 7/3 — left panels — and j/4 — right panels — eigenfunctions for the (/-modes of a typical 0.6 Mq 
white dwarf — top panels — for BPM 37093 — central panels — and for PG 1159-035 — bottom panels — as a 
function of the mass coordinate log(l — m/M*). For the case of BPM 37093 the crystallized core is shown as a hatched 
area. See text for details. 



where R± and Ro are the stellar radius and the radial co- 
ordinate of the crystallization front, respectively. By using 
Poisson's equation 



r 2 dr V dr 



(5) 



and Eq. (4) it can be shown after a straightforward calcu- 
lation that the dimensionless strain, hj^f , is given by 



n 33 



5 x 10 



-18 



R* \ 2 

Rq) v ■ i inllz 



A{R+) - F^F 2 R A{R )\ cos(at) 
being v — 2ira the frequency of the signal and 



(6) 



A{r) 

2/3 



2/4 - 2y 3 
gr 

l d& 
9 dr 



(7) 



F R = ^ 



Mo 

Ro 
R~* 



In these expressions $' is the perturbed gravitational po- 
tential, g is the gravitational acceleration at a given ra- 
dius, M* and i?* are the mass and radius of the star, and 
M and Ro are the mass and radius of the crystallized 
core of the white dwarf. Obviously, in the non-crystallized 
case, Mo and Rq are identically zero (corresponding to the 
stellar center). The advantage of the previously described 
formalism is that the quantities 1/3 and j/4 can be easily 
tabulated for a typical stellar model, whereas the rest can 
be observationally obtained. Finally, the luminosity radi- 
ated in the form of gravitational waves is (Osaki & Hansen 
1973): 



10 



M Q 



! (—Y(- 

\RqJ vi 



mHz 



A(R+) - F^AiRo) 



(8) 
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Fig. 2. Same as Fig. 1 for the /-mode of our 0.6 Mq fiducial model, for BPM 37093 and for PG 1159-035. 



5. Results 

Figure 1 shows the run, as a function of the mass co- 
ordinate, of the functions ys — left panels — and 2/4 - 
right panels — discussed in §4 for the g-modes of our fidu- 
cial model (a typical 0.6 Mq white dwarf made of carbon 
and oxygen with a liquid core and hydrogen envelope of 
10- 4 M„), for BPM 37093 — a massive (M* ~ 1.05 Mq) 
white dwarf, with a sizeable crystallized core — and for 
PG 1159-035, the other known white dwarf with unam- 
biguously identified quadrupole (/-modes. It is important 
to mention at this point that the fractional change in ra- 
dius due to pulsations, SR^/R^, must not necessarily be 
the same for each pulsation mode. However, the lineal 
theory of nonradial pulsations does not provide any clue 
about the value of the fractional change in radius, since the 
governing equations are homogeneous and the normaliza- 
tion of the eigenfunctions is arbitrary (Cox 1980). In ad- 
dition, 5R+/R+ is poorly constrained by the observations, 
since the luminosity variations of pulsating white dwarfs 
are almost exclusively caused by changes in temperature, 
not by their radius variations (Robinson et al. 1982). Thus, 
we have adopted, somehow arbitrarily, that the fractional 
change in radius due to pulsations is 5R+/R+ = 10~ 4 for 
all the considered modes, which is a typical value for pul- 



sating white dwarfs, and reproduces reasonably the am- 
plitude of the observed light curve (Robinson et al. 1982). 

We start discussing the results obtained for our fidu- 
cial model (top panels of Fig. 1). We have chosen to dis- 
play the quadrupole (£ = 2) <?-mode with radial order 
k = 25, which has a period P = 678.22 s. As can be 
seen the functions j/3 and j/4 are small everywhere in the 
star, being their amplitudes of the order of ~ 10~ 6 and 
~ 10~ 4 , respectively. Moreover, the amplitudes are only 
significant for the central regions of the white dwarf. For 
the case of BPM 37093 — central panels — we show the 
quadrupole g-mode with k = 27. This g-mode has a pe- 
riod P = 536.4 s, which is very close to one of the observed 
periods, which has P = 531.1 s. It is important to realize 
that since BPM 37093 is a massive white dwarf, a size- 
able region of its core is crystallized. This region is clearly 
marked in the central panels of Fig. 1 as a shaded area. 
We note that in this region the amplitudes of both j/3 and 
i/4 are null. Finally, for the PG 1159-035 model — bottom 
panels — we display the quadrupole mode with k = 30, 
which best fits the observed period of P = 423.2 s. This 
mode has a period P = 423.8 s, thus providing an excel- 
lent fit to the observational data. 
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Fig. 3. Same as Fig. 1 for the p-modes of our 0.6 M Q fiducial model, for BPM 37093 and for PG 1159-035 

Table 1. Summary of the gravitational wave emission of the g-modes of BPM 37093 and PG 1159-035. Our fiducial 
model is also shown for the sake of comparison. In all cases we have adopted 6Rt/R+ = 10~ 4 . 



Model 


M/Mq 


k 


Po (s) 


P c (s) 






Law (erg/s) 


log(E K ) (erg) 


BPM 37093 


1.10 


26 


511.7 


516.7 


4.8 x 10" 


28 


9.2 x 10 18 


46.4 






27 


531.1 


536.4 


5.4 x 10" 


28 


1.1 x 10 19 


46.5 






28 


548.4 


555.8 


5.3 x 10" 


28 


9.9 x 10 18 


46.6 






29 


582.0 


574.9 


6.4 x 10" 


28 


1.3 x 10 19 


46.7 






30 


600.7 


593.0 


6.4 x 10" 


28 


1.2 x 10 19 


46.8 






32 


633.5 


630.4 


5.5 x 10" 


28 


8.3 x 10 18 


46.9 


PG 1159-035 


0.54 


25 


352.7 


358.9 


1.0 x 10" 


30 


5.0 x 10 Ib 


44.0 






30 


423.8 


423.2 


2.3 x 10" 


30 


1.7 x 10 17 


43.8 






50 


694.9 


684.5 


3.5 x 10" 


31 


1.6 x 10 15 


43.4 






55 


734.2 


752.9 


2.0 x 10" 


31 


4.1 x 10 14 


43.3 






60 


812.6 


818.1 


1.8 x 10" 


31 


2.8 x 10 14 


43.2 






70 


968.7 


950.1 


6.0 x 10" 


32 


2.5 x 10 13 


42.8 


0.6 M Q 


0.6 


1 




66.6 


6.9 x 10" 


25 


1.0 x 10 28 


47.0 




10 




310.3 


1.8 x 10" 


27 


3.1 x 10 21 


44.8 






20 




555.2 


1.5 x 10" 


27 


6.6 x 10 20 


45.1 



In Table 1 we summarize the most important results 
for several (/-modes of the models computed so far. The 
first column of this table lists the model. The second col- 
umn corresponds to its respective mass. In the third col- 
umn we show the radial order, k, of the computed g-mode. 



The observed and the computed periods (in seconds) are 
given in columns 4 and 5, respectively. We have assumed 
that all of the observed periods in BPM 37093 and in 
PG 1159-035 in Table 1 are I = 2, following the works 
by Kannan et al. (2005) and Winget et al. (1991), re- 
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spectively. In column 6 we give the maximum dimension- 
less strain, h max , as computed from Eqs. 6 and 7. For 
BPM 37093 we have computed it adopting the measured 
distance to the source (d = 16.8 pc), whereas for our fidu- 
cial model we have adopted a distance d — 50 pc, which 
we consider to be representative of a typical white dwarf. 
For the case of PG 1159-035 we have adopted a distance 
of 400 pc, in line with the determinations of Werner et al. 
(1991) and Kawaler & Bradley (1994). Column 8 provides 
the luminosity in the form of gravitational waves radiated 
away, Low, computed with Eq. (8). Finally, in the last 
column of Table 1 we list the kinetic energy of each one of 
the modes. Note that, in general, the agreement between 
the computed and the observed periods is rather good for 
all the modes, both for the case of BPM 37093 and for 
PG 1159-035. However, the amplitudes of the dimension- 
less strains are rather small in all the cases. This is also 
the case for the luminosities radiated away in the form of 
gravitational waves. In the best of the cases BPM 37093 
radiates away <~ 10 19 erg/s in the form of gravitational 
waves, whereas PG 1159-035 radiates away a much more 
modest amount, only ~ 10 17 crg/s. 

It is interesting to note that for our fiducial model — 
third section of Table 1 — the larger the radial order k, 
the smaller the dimcnsionless strain is and the smaller are 
the luminosities radiated away in the form of gravitational 
waves. In particular, an increase from k = 1 to k = 10 
produces a reduction of a factor of almost 4 x 10 2 in the 
dimensionless strain and of 3 x 10 6 in the flux of gravita- 
tional waves. The reductions when considering the k = 20 
mode are much more modest. The reason for this is the 
following, low-fc modes sample the core more than high-Zc 
modes, and since the core has a higher density, larger mass 
motions are produced, and hence more gravitational wave 
losses are produced. Nevertheless, both the dimcnsionless 
strains and the fluxes of gravitational waves are in this 
case much larger than those found for BPM 37093 and 
PG 1159-035. The reasons for this behavior are easy to 
understand. For the case of BPM 37093 pulsations occur 
only in a small region of the star as a result of its crystal- 
lized core. Thus, despite its mass being much larger than 
that of our fiducial model the emission of gravitational 
waves is strongly inhibited. For the case of PG 1159-035, 
the most important reason why so few gravitational waves 
are radiated away is its small mass (and average density). 

Given the results obtained for the quadrupole g-modes 
studied up to now it is natural to wonder whether other 
modes, namely the /- and p-modes, of pulsating white 
dwarfs can radiate away a measurable amount of gravita- 
tional waves. To this regard we have extended our calcu- 
lations to incorporate such modes, despite the lack of ob- 
servational evidence for them. Obviously, for these modes 
we do not know which is the appropriate value of 5R+/ R+, 
since the estimate SR+/R* = 10~ 4 is based on observed 
g-modes in white dwarfs. However, for the calculations re- 
ported here we adopted the same value. Fig. 2 shows the 
run of the functions and 1/4 for /-modes of our fiducial 
model, for BPM 37093 and for PG 1159-035. As is the 



0.6 M 




10- 5 10- 4 10- 3 io- s 10- 1 10° 10 1 10 2 



i/(Hz) 

Fig. 4. A comparison of the signal produced by the 
quadrupole g-modes of our fiducial model — circles - 
by the /-mode — square — and by the p-modes — tri- 
angles — with the spectral distribution of noise of LISA 
for a one-year integration period, and assuming that the 
source is located at 50 pc. 

case for all /-modes the radial order is zero, and their re- 
spective periods are (from top to bottom): 11.35 s, 3.7 s 
and 20.9 s, respectively. We note that now the functions 2/3 
and 2/4 are much larger than in the case previously studied. 
Moreover, it is worth noting that in this case the functions 
2/3 and 2/4 have large amplitudes everywhere and do not 
vanish at the surface. Consequently, we expect that a large 
amount of gravitational waves can be radiated away. This 
is indeed the case, as can be observed in Table 2, where 
we show the model (first column), its respective mass (sec- 
ond column), the considered mode (third column), the ra- 
dial order (fourth column), the respective frequency of the 
mode (fifth column), the dimcnsionless strain (sixth col- 
umn), the luminosity radiated away in the form of gravita- 
tional waves (seventh column) and the kinetic energy (last 
column). Note that for all the /-modes of the three models 
presented here the dimensionless strains are several orders 
of magnitude larger and, moreover, the luminosities radi- 
ated away are much larger than the optical luminosities, 
even of the order of 10 41 erg/s in the case of BPM 37093. 
This, in turn, could have important consequences since it 
could provide one of the possible reasons why these modes 
have not been observed thus far: if they are excited they 
are quickly damped by emission of gravitational waves. 

Now we turn our attention to the p-modes. Again, we 
show the run of the 2/3 and 2/4 functions of the p-modes 
in terms of the mass coordinate for our three models in 
Fig. 3. We have chosen to show the k = 25 mode for 
all three cases. Their respective periods are P = 0.56 s for 
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Table 2. Summary of the gravitational wave emission of the /- and p-modes of BPM 37093 and PG 1159-035. Our 
fiducial model is also shown for the sake of comparison. In all cases we have adopted 5R+/R+ = 10~ 4 . 



Model 


M/M q 


Mode 


k 


v (Hz) 








L G w (erg/s) 


\og{E K ) (erg) 


BPM 37093 


1.10 


/ 





2.7 x 10" 


i 


6.6 x 10" 


19 


3.4 x 10 4i 


49.2 






P 


1 


9.3 x 10" 


i 


9.1 x 10" 


20 


7.7 x 10 40 


47.8 


PG 1159-035 


0.54 


f 





4.8 x 10" 


•2 


1.2 x 10" 


Tl 


2.1 x 10 3b 


45.5 






P 


1 


5.9 x 10" 


2 


7.1 x 10" 


23 


5.1 x 10 34 


44.6 


0.6 M 


0.60 


f 





8.8 x 10" 


! 1 


5.6 x 10" 


20 


2.3 x 10 3a 


47.0 




P 


1 


1.8 x 10" 


1 


6.4 x 10" 


21 


1.3 x 10 38 


48.6 






P 


5 


4.6 x 10" 


1 


6.1 x 10" 


21 


7.5 x 10 38 


46.0 



our fiducial 0.6 M Q model, P = 0.11 s for BPM 37093 and 
P = 2.26 s for PG 1159-035. We note that the amplitudes 
of the y3 and j/4 functions for p-modes are much smaller 
than those of the corresponding /-modes, and comparable 
to the corresponding g-modes studied before. In addition, 
in contrast with the situation for the /-modes, the ampli- 
tudes of 2/3 and 2/4 are almost negligible in regions close to 
the surface of the white dwarf. However, because the pul- 
sation frequencies of p-modes are considerably higher than 
those of the /- and g-modes, the dimensionless strains 
- see Table 2 — are consequently large and the corre- 
sponding gravitational wave luminosities are very large 
as well, although roughly one order of magnitude smaller 
than those obtained for the /-modes. 

In order to check whether or not LISA would be able 
to detect the pulsating white dwarfs studied here we have 
proceeded as follows. We first have assumed that the inte- 
gration time of LISA will be one year. The signal-to-noise 
ratio, 77, is given by: 



BPM 37093 



v 



h 2 (a) da 
S(a) 2^ 



(9) 



where S(cr) = Sh(a)T is the sensitivity of LISA, r is 
the integration period, h(a) is the Fourier Transform 
of the dimensionless strain, and a has been previously 

defined. It can easily be shown that for a monochromatic 

1/2 

gravitational wave r\ = h(a)/S h (a). We have adopted 
a signal-to-noise ratio r\ = 5. We have furthermore used 
the integrated sensitivity of LISA as obtained from 
http : //www . srl . caltech . edu/^shane/sensitivity 
The results of this procedure are shown in Figs. ^ 
and for our fiducial model (0.6 Mq carbon-oxygen 
white dwarf), for BPM 37093 and for PG 1159-035, 
respectively. In all three figures the g-modes are shown 
as circles, the /-mode is shown as a square and the 
p-modes are displayed as triangles. As previously done, 
for our fiducial model we have adopted a distance of 
50 pc, for PG 1159-035 we have assumed a distance of 
400 pc, and for BPM 37093 we have used its measured 
distance (16.8 pc). As can be seen, for none of the three 
cases studied here LISA will be able to measure the 
dimensionless strains of the g-modes, even at a reduced 
signal-to-noise ratio. Most of the p-modes will not be 
observed as well, either because their frequencies are too 
high to be observed by LISA or because they are too 




10- 5 10- 4 10" 3 10" 2 10- 
i/(Hz) 



Fig. 5. Same as Fig. 4 for the case of BPM 37093. The 
distance in this case is known, d = 16.8 pc. 



weak. Moreover, given that these modes radiate huge 
amounts of energy in the form of gravitational waves they 
become quickly damped and, consequently, it will be very 
difficult to detect them. Particularly, and given that for 
the p-modes studied here the time-averaged dissipation 
rates of pulsations due to radiative (photons) heat 
leakage and neutrino losses are much smaller than the 
luminosity radiated as gravitational waves, an estimate 
of the damping timescale, r^, can be easily obtained by 
considering the kinetic energy of the mode, E^, which 
can be easily computed from our numerical models: 



2E 



Td 



K 



L 



GW 



(10) 



For instance, for the p-mode with k = 1 of BPM 37093 
one obtains Td — 0.5 yr, which clearly is too short to allow 
a detection. 

Finally, the three /-modes of the models presented 
here lay in the appropriate range of frequencies and, addi- 
tionally, they are well over the sensitivity curve of LISA. 
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Fig. 6. Same as Fig. 4 for the case of PG 1159-035. Since 
the distance to PG 1159-035 is not accurately known we 
have adopted d = 400 pc, which is a reasonable estimate. 

However, as it was the case for the p-modes, they also 
radiate very large amounts of gravitational waves and, 
hence, they will be quickly damped, hampering the pos- 
sibility of detection. In this case we obtain a damping 
timescale r d ~ 2.9 yr for the /-mode of BPM 37093. Note 
that although the luminosity radiated away in the form 
of gravitational waves is larger for the /-mode than for 
the p-mode considered previously the damping timescale 
is larger. This is so because the kinetic energies involved 
are quite different: Ek — 1-54 x 10 49 erg for the /-mode 
and Ek = 6.42 x 10 47 erg for the p-mode with k = 1. 
For the sake of completeness we present in Table 2 all the 
relevant information for the p- and /-modes which could 
be eventually detected. 

6. Discussion and conclusions 

In this paper we have computed the gravitational wave 
emission of pulsating white dwarfs. We have started 
by computing the gravitational wave radiation of white 
dwarfs undergoing nonradial g-raode pulsations, which are 
currently observed in a handful of pulsating white dwarfs. 
We have focused on three model stars. Our fiducial model 
corresponded to an otherwise typical 0.6 Mq model white 
dwarf with a carbon-oxygen fluid core and a hydrogen 
envelope. We have also paid attention to two additional 
white dwarf models, corresponding to two stars for which 
quadrupole g-modes have been observed so far, namely, 
BPM 37093 and PG 1159-035. We have shown that in 
all these cases the gravitational wave signal is too weak 
to be observed by future space-borne interferometers, like 
LISA. More specifically, we have found that the luminosi- 



ties in the form of gravitational waves radiated away by 
these stars and the corresponding dimensionless strains 
are very small in all the cases, in agreement with the pio- 
neering work of Osaki & Hansen (1973). Hence, all these 
sources contribute to the Galactic noise and no individ- 
ual detection are expected, despite the proximity of the 
sources. For the sake of completeness we have also elabo- 
rated on this subject and we have computed the gravita- 
tional wave emission of white dwarfs undergoing nonradial 
/- and p-mode oscillations, even if these modes have not 
been observationally detected whatsoever. We have found 
that for white dwarfs undergoing these kind of pulsations 
the luminosities in the form of gravitational waves radi- 
ated away are very large in all the cases, in line with the 
earlier results of Osaki & Hansen (1973). Consequently, 
these modes, if excited, should be very short-lived, thus 
hampering their eventual detection. 

It may seem that for the case of pulsating white dwarfs 
undergoing g-mode oscillations there could still be a pos- 
sibility of indirect detection by measuring the secular rate 
of change of the period of the observed modes. However, 
this is not the case. Particularly, it can be easily shown 
that the secular rate of change of the period of a pulsating 
white dwarf is given by 



P T R+ 

p = - a T +b R: 



(11) 



where T is the temperature of the isothermal core and 
a and b are constants of order unity which depend on 
the chemical composition, thicknesses of the H atmosphere 
and He buffer, equation of state, and other ingredients 
involved in the modeling of white dwarfs. For DA white 
dwarfs in the ZZ Ceti instability strip, the second term 
of the right hand side of Eq. (|llf> is usually negligible 
and, thus, the secular rate of change of the period only 
reflects the speed of cooling — see, for instance, Isern et 
al. (1992) and references therein — whereas for GW Vir 
stars this term is relevant, but can be accounted for by the 
theoretical models and, hence, the speed of cooling can 
also be derived. However, any additional source of extra 
cooling — like gravitational waves — would eventually 
translate into an anomalous rate of period change (Isern 
et al. 1992): 



P o . _ Lqw 
P c ~ L + L v 



(12) 



where P is the observed period, P c is the computed pe- 
riod without taking into account the emission of gravi- 
tational waves, L v is the neutrino luminosity — which 
is important for hot white dwarfs — and the rest of the 
symbols have been previously defined. The rate of sec- 
ular period change has been measured for some pulsat- 
ing white dwarfs. Particularly, for G117-B15A (Kepler 
et al. 2000) it has been possible to measure the sec- 
ular variation of the main observed period of 215.2 s, 
P = (2.3 ± 1.4) x 10~ 15 s s -1 , with unprecedented ac- 
curacy. In fact, this white dwarf is the most stable optical 
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clock known, and has been used to pose tight constraints 
on the mass of the axion (Corsico et al. 2001b) and the 
rate of variation of gravitational constant (Benvenuto et 
al. 2004). Other pulsating white dwarfs — like L 19-2 and 
R 548 — have also determinations of the secular rate of pe- 
riod change but are not as accurate as that of G117-B15A. 
Note, however, that for a 0.6 Mq white dwarf undergoing 
quadrupole g-mode oscillations with a period of P ~ 200 s 
the right hand side of Eq. I)12|l amounts to 2 x 10 -10 and 
P c ~ 10~ 15 s s _1 . In other words the gravitational radi- 
ation will produce a change in P ~ 10~ 25 sr 1 , making 
impossible such an indirect detection even if accurate ob- 
servational data and reliable theoretical models become 
eventually available. Again, the only case of interest here 
would be the case in which Lqw ~ L + L u , which may be 
true for the /- and p-mode pulsators. 
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